Endohedral Impurities in Carbon Nanotubes 
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A generalization of the Anderson model that includes pseudo-Jahn- Teller impurity coupling is 
proposed to describe distortions of an endohedral impurity in a carbon nanotube. Treating the 
distortion within mean-field theory, spontaneous axial symmetry breaking is found when the vibronic 
coupling strength g exceeds a critical value gc. The effective potential in the symmetry-broken state 
is found to have 0(2) symmetry, in agreement with numerical calculations. For metallic zigzag 
nanotubes endohedrally-doped with transition metals in the dilute limit, the low-energy properties 
of the system may display two-channel Kondo behavior; however, strong vibronic coupling is seen 
to exponentially suppress the Kondo energy scale. 

PACS numbers: 61.46.-l-w, 64.70.Nd, 73.22.-f, 73.63.Fg 



Single wall carbon nanotubes (SWNTs) are often 
cited as having the potential to contribute to a num- 
ber of important technologies, such as electrochemical 
storage, nanoelcctronics, medical imaging, and quantum 
computing. In these examples, it is the controlled intro- 
duction of impurities, either exohedrally or endohedrally, 
that are of importance. Experimentally, as with the 
fuUerenes, SWNTs have been doped by a variety of atoms 
and molecules. Recent techniques have led to "peapod 
nanotubes," where one dimensional arrays of fuUerenes 
P and metallofuUerenes |^ have been encapsulated 
by SWNTs. It was subsequently demonstrated Q using 
scanning tunneling microscopy (STM) that the endohe- 
dral doping of SWNT with fuUerenes does indeed modify 
the local electronic structure of the SWNT. 

For the case of SWNTs endohedrally-doped with atoms 
or small molecules. Cole and collaborators |^ have identi- 
fied two distinct structural phases: the axial phase, where 
the impurity is confined to the vicinity of the tube's axis 
of symmetry, and the shell phase, where the impurity 
is confined to an internal cylindrical surface concentric 
with the tube. For Uthium in a (5,5) SWNT, the stabil- 
ity of the shell phase has been seen numerically , while 
the axial phase is the ground state for a system of II2 
molecules at low densities in a (5,5) SWNT [Q. 

Here, the focus is on the development of a theoretical 
framework for describing a single impurity in a SWNT. 
The shell phase can be thought of a resulting from the 
spontaneous breaking of axial symmetry within mean- 
field theory. Such a description is the ferroelastic analog 
of local magnetic moment formation in a metal within 
the Anderson model. The inversion symmetry breaking 
in metallofuUerenes proceeds along similar lines ||]. 

It is anticipated that, in analogy with magnetic im- 
purities in non-magnetic metallic hosts, the dynamics of 
the vibronic impurity in a metallic nanotube will lead 
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FIG. 1: Sketch of a section of a (5, 5) SWNT with an endo- 
hedral atom off-axis (shell phase). 



to Kondo-related effects. In a different context, spin-flip 
Kondo effects in a carbon nanotube system already have 
been experimentally explored [ pO| where the electronic 
states of a section of a nanotube, analogous to the local 
electronic states of a quantum dot, play the role of the 
impurity states. Kondo effects for magnetic Co clusters 
adsorbed on metallic SWNT have been discovered 
using low temperature STM and have been subsequently 
analyzed theoretically fl^ . 

The vibronic states resulting from endohedral doping 
and the subsequent Kondo coupling resulting from vir- 
tual occupation of those levels give rise to another kind 
of possible Kondo system. For strong vibronic coupling, 
the Kondo temperature, roughly the temperature below 
which Kondo-related effects are manifested, is exponen- 
tially suppressed. 

There are three classes of electronic band structure for 
the SWNT and they can be characterized by the 
degeneracy at the Fermi energy. The first class is semi- 
conducting with no density of states at the Fermi energy. 
The remaining SWNTs are metallic with the second class 
having two non-degenerate bands crossing at the Fermi 
energy away from the zone center; the third class has two 
doubly degenerate bands crossing at the zone center at 
the Fermi energy. The armchair SWNTs (n, n) fall in the 
second class, and the zigzag SWNTs {3q, 0) q G fall 
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in the third class, while all others are semiconducting. 

For the case where the SWNT has a single non- 
degenerate band at or near the Fermi energy (the sec- 
ond class and many cases of the first class) , the aig band 
electrons may hop to the localized electronic state of the 
same parity on the impurity site. This impurity state 
is taken to be vibronically coupled through the linear 
pseudo-Jahn- Teller effect to a low-lying state of oppo- 
site parity under the transverse displacement Q of the 
impurity from the axially symmetric position. (Actually, 
this displacement can be regarded as a normal coordinate 
that includes an accompanying distortion of a nanotube 
unit cell with the appropriate e2„ symmetry.) Vibronic 
terms that are of 0{Q^) and higher are neglected here, as 
they have been seen to be small in the related endohedral 
fuUcrene systems |^. 

Thus, the electronic Hamiltonian is taken as 



■?ie = ^ ^kclck + ^ Emd]^dm + ^ (Vfecj,do + H 



(1) 

where m — Q, ±1 corresponding to the axial angular mo- 
mentum of the electronic states and the localized one- 
electron energies have a double degeneracy with Ei = 
E^i = E^. The operators c\ (cfe) create (destroy) an 
electron on the nanotube with wavevector k (fc > 0). 
The operator dj^ {dm) creates (destroys) an electron on 
the unit cell containing the endohedral dopant. All en- 
ergies are relative to the Fermi energy. The Hamiltonian 
is diagonal in the spin index. Thus, spin indices are sup- 
pressed. 

The discussion here is restricted to the case where the 
on-site Coulomb repulsion U is sufficiently large that in 
the limit Vk — Q only the singly occupied m = impurity 
state is below the Fermi energy Ef (infinite-C/ model in 
the local moment regime). Estimates based on metallo- 
fuUcrene calculations indicate this case to be relevant 
to alkali-doped nanotubes. 

The pseudo-Jahn- Teller contribution to the Hamilto- 
nian must be invariant under symmetry operations of 
the nanotube. For a nanotube with at least D3 symme- 
try, the vibronic term for this model has continuous axial 
symmetry. The form of the coupling term is determined 
by symmetry up to an overall coupling constant g to 
be 



WjT = -5 X! {Q-md\ndo + H.C.) 



(2) 



m=±l 



where Q±i ^ Qx ± iQy 

Lastly, the form of the elastic energy of the distortion 
is taken to be 



(3) 



m— ±1 



This model is now considered for a classical field Q, 
a mean-field theory in the distortion. Writing the dis- 
tortion in polar form Qm = QeKp{im(f)), the phase is 



gauged away by transforming dm —>■ dm exp(— im0), leav- 
ing only the magnitude Q. The Hamiltonian is bilinear in 
operators and thus is easily diagonalized. The impurity 
Green's function Gqo is simply that of the non-interacting 
Anderson model with an additional contribution to the 
self-energy due to the pseudo-Jahn- Teller term 



Goo = (0|G|0) = 



1 
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with the self-energy S(a;) 



' + ir] - Ek LO - E^ 



(4) 



(5) 



The energy change due to the vibronic impurity is 
given by 



E=- [ dwImlnGoo(w) + 
n J 2 



(6) 



For systems with coupling strength g > gc, a distorted 
equilibrium with Q ^ is energetically favorable. All 
distortion phases 4>, however, are degenerate, correspond- 
ing to an adiabatic potential with an equipotential circle 
concentric with the symmetry axis. Such 0{2) invariance 
is broken with the inclusion of higher-order Jahn- Teller 
couplings. As numerical calculations |^ point to con- 
tinuous rotational symmetry of the equilibrium positions 
of the dopants considered, the linear PJT model proposed 
is consistent with these results. 

In analogy with the Anderson model for magnetic im- 
purities, the broken symmetry is an artifact of the mean- 
field description. Phase fluctuations of the distortion 
stemming from zero point motion of the dopant restore 
axial symmetry. 

For a linear conduction band with width 2t, Vk inde- 
pendent of k and Eq near Ej, it is straightforward to 
derive the Stoner-like criterion that 
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and 



/(7,e*,eo) = / d£_— 

Jo ? + e* 
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(7) 



(8) 



where 7 is the width of the virtual bound state resonance 
near Ef in units of t, 7 = 7r|ypp/i, e* = and 
cq — {Eq — e)/t where e is the shift in the resonance at 
Ef due to the hybridization self-energy. The integral in 
Eq. H can be expressed in terms of elementary functions, 
yielding the analytic form for Eq. |^ which is plotted in 
Fig. I 

For the case where the SWNT has a doubly degenerate 
band at or near the Fermi energy (the third class and the 
remainder of the first class) , d-states of a transition metal 
impurity could readily hybridize with the Cg bands of the 
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FIG. 2: Phase diagram of vibronic impurity in a SWNT. is 
measured in units of fct, 7 is the resonance width in units of 
t. Phase boundaries are plotted for e* = 0.2, and eo = —0.05 
(sohd curve), eo = —0.1 (dashed), eo = —0.15 (thick), and 
eo = -0.25 (dotted). 

nanotube, leading to an electronic model of the following 
form 

+ Ekm{VkcLd„,,+R.c.) (9) 

where p = ±1 is the parity of the localized orbital dou- 
blets. 

The impurity states of opposite parity can be coupled 
through the impurity distortion via the pseudo-Jahn- 
Teller effect. 

'HjT = -5 51 {rnQrndl^dm-i + H.c.) (10) 

m 

As in the case of the Hamiltonian in Eq. ||, the phase of 
the displacement can be gauged away with the transfor- 
mation 

dmp ^ d^pe^S^^^+t™-!)^) 

The pseudo-Jahn- Teller contribution to the self-energy in 
this {Eu (B Eg) ^ Cu model is found to be one-half that 
of Eq. Thus, a qualitatively similar phase diagram is 
found for the two-band model with g"^ ^ g'^/2. 

For SWNTs of the third class, for g < gc, orbitally 
degenerate conduction electrons are coupled to orbital 
impurity states through virtual fluctuations in impurity 
occupancy. Orbital states can be labeled by a pseu- 
dospin, leading to a possible two-channel Kondo descrip- 
tion where the real spin gives rise to the channel indices 
p6[ ; however, the low density of states p of metallic 
zigzag SWNTs tends to lower the Kondo temperature 
Tk relative to standard metallic hosts. 



Under linear vibronic mixing, the ground state double 
degeneracy of the isolated impurity is preserved. For fi- 
nite U, provided that fluctuations between the ground 
state orbital doublet and a doubly occupied, doubly de- 
generate set of states dominate over fluctuations to the 
unoccupied singlet, the two-channel Kondo description is 
still operative; however, Tk is modified by changes in the 
hybridization and impurity energy levels. 

This can be quantified in the strong coupling limit 
g 3> gc- At low temperatures, fluctuations of the mag- 
nitude of the distortion can safely be neglected. Hence, 
the impurity atom acts as a quantum rotor, displaced 
from the tube axis by Qo, the mean-field value of the 
magnitude of the distortion. 

The hybridization term can be rewritten as 

Ha-5^(T4cL|0)(m|+H.c.) (12) 

km 

where |to), (m = ±1) are the impurity eigenstates with 
energy ^ = Ei + ^{6 ~ A), S ^ E^i ~ Ei and = 
S'^+4g^Qg. The hybridization is reduced by vibronic mix- 
ing in the strong coupling regime to give Vk = Vk ^^^^ ■ 

Additionally, with increasing coupling strength, the vi- 

2 

bronic energy levels drop, with ^ ^ — ^ as 5 — > 00. 

Collectively, this implies that the Kondo exchange 
coupling J, as obtained from the generalized Ander- 
son model through a Schrieffer- Wolff transformation, be- 
haves asymptotically as J ~ as g — > 00, leading 
to Txig) ^ TxiO) exp{-g'^ /V^Kp). Strong vibronic cou- 
pling can quickly suppress the temperature scale below 
which Kondo effects will be manifested. 

In summary, a vibronic generalization of the Anderson 
model is seen to be a convenient framework for under- 
standing the structural phases of endohedral impurities 
in SWNTs. An analogy is noted between the formation of 
impurity moments in the magnetic Anderson model and 
the distortion of the endohedral dopant in this vibronic 
generalization. 

Within mean-field theory, the high symmetry axial 
phase exists for vibronic coupling g < gc- For strong 
coupling, the shell phase is seen to be stable. Beyond 
mean-field theory, for the models considered where or- 
bital degeneracy becomes vibronic degeneracy, a modifi- 
cation of the two-channel Kondo effect is operative; how- 
ever, Tk is seen to be exponentially suppressed for strong 
vibronic coupling. 
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